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Abstract 

In papers[l, 2], the amplitudes of continuous-time quantum walk on graphs possessing 
quantum decomposition (QD graphs) have been calculated by a new method based on 
spectral distribution associated to their adjacency matrix. Here in this paper, it is shown 
that the continuous-time quantum walk on any arbitrary graph can be investigated by 
spectral distribution method, simply by using Krylov subspace-Lanczos algorithm to 
generate orthonormal bases of Hilbert space of quantum walk isomorphic to orthogonal 
polynomials. Also new type of graphs possessing generalized quantum decomposition 
have been introduced, where this is achieved simply by relaxing some of the constrains 
imposed on QD graphs and it is shown that both in QD and GQD graphs, the unit 
vectors of strata are identical with the orthonormal basis produced by Lanczos algorithm. 
Moreover, it is shown that probability amplitude of observing walk at a given vertex is 
proportional to its coefficient in the corresponding unit vector of its stratum, and it can be 
written in terms of the amplitude of its stratum. Finally the capability of Lanczos-based 
algorithm for evaluation of walk on arbitrary graphs ( GQD or non-QD types), has been 
tested by calculating the probability amplitudes of quantum walk on some interesting 
finite (infinite) graph of GQD type and finite (infinite) path graph of non-GQD type, 
where the asymptotic behavior of the probability amplitudes at infinite limit of number 
of vertices, are in agreement with those of central limit theorem of Ref.[3]. 

Keywords: Continuous-time quantum walk, Spectral distribution, Graph, 
Krylov subspace, Lanczos algorithm. 

PACs Index: 03.65.Ud 
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Random walks on graphs are the basis of a number of classical algorithms. Examples include 
2-SAT (satisfiability for certain types of Boolean formulas), graph connectivity, and finding 
satisfying assignments for Boolean formulas. It is this success of random walks that motivated 
the study of their quantum analogs in order to explore whether they might extend the set of 
quantum algorithms. 

Recently, the quantum analogue of classical random walks has been studied in a flurry of 
works [4, 5, 6, 7, 8, 9]. The works of Moore and Russell [8] and Kempe [9] showed faster 
bounds on instantaneous mixing and hitting times for discrete and continuous quantum walks 
on a hypercube (compared to the classical walk). 

A study of quantum walks on simple graph is well known in physics(see [10]). Recent studies 
of quantum walks on more general graphs were described in [4, 5, 7, 11, 12]. Some of these 
works studies the problem in the important context of algorithmic problems on graphs and 
suggests that quantum walks is a promising algorithmic technique for designing future quantum 
algorithms. One approach for investigation of continuous-time quantum walk on graphs is 
using the spectral distribution associated with the adjacency matrix of graphs. Authors in 
[1, 2] have introduced a new method for calculating the probability amplitudes of quantum 
walk based on spectral distribution, where a canonical relation between the Hilbert space of 
stratification corresponding to the graph and a system of orthogonal polynomials has been 
established, which leads to the notion of quantum decomposition (QD) introduced in [13, 14] 
for the adjacency matrix of graph. Also it is shown in [1] that by using spectral distribution 
one can approximate long time behavior of continuous-time quantum walk on infinite graphs 
with finite ones and vice versa. In [1, 2], only the particular graphs of QD type have been 
studied . 

Here in this work we try to investigate continuous-time quantum walk on arbitrary graphs 
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by spectral distribution method. To this aim, first by turning the graphs into a metric space 
based on distance function, we have been able to generalize the stratification and quantum 
decomposition introduced in [13], such that the basis of Hilbert space of quantum walk consist 
of superposition of quantum kets of vertices belonging to the same stratum, but with different 
coefficients, while the coefficients are the same in QD case, therefore QD graphs introduced 
in [1, 13] are particular kind of graphs possessing generalized quantum decomposition (GQD). 
Then we show that both in QD and GQD graphs, the unit vectors of strata are identical 
with the orthonormal basis produced by Lanczos algorithm. Also, in the case of GQD graphs 
we show that probability amplitude of observing walk at a given vertex is proportional to 
its coefficient in corresponding unit vector of its stratum, and it can be written in terms of 
the amplitude of its stratum. For more general graphs, the Lanczos algorithm transforms the 
adjacency matrix into a tridiagonal form (quantum decomposition) iteratively, where we use 
this fact for studying non-QD type graphs. Indeed, the Lanczos algorithm gives a three-term 
recursion structure to the graph, so the spectral distribution associated to adjacency matrix 
can be determined by Stieltjes transform. In order to see the power of Lanczos-based algorithm 
in the investigation of continuous-time quantum walk on arbitrary graphs ( GQD or non-QD 
types), we have calculated the amplitudes of quantum walk on some interesting finite (infinite) 
graph of GQD type and finite (infinite) path graph of non-GQD type. 

The organization of the paper is as follows. In section 2, we review the Krylov subspace 
methods and Lanczos algorithm. In Section 3, we give a brief outline of some of the main 
features of graphs and introduce generalized stratification. Section 4 is concerned with the 
Hilbert space of generalized stratification. In Section 5, we review the Stieltjes transform 
method for obtaining spectral distribution /i, and establish an isometry between orthogonal 
polynomials and Hilbert space of generalized stratification. Section 6 is devoted to the method 
for computing amplitudes of continuous-time quantum walk, through spectral distribution \i 
of the adjacency matrix A. In section 7 we calculate the amplitudes of quantum walk on some 



continuous-time Quantum walk 



5 



interesting finite (infinite) graph of GQD type and finite (infinite) path graph of non-GQD 
type. At the end we study the asymptotic behavior the probability amplitudes at infinite limit 
of number of vertices, where the results thus obtained are in agreement with those of central 
limit theorem of Ref.[3]. Paper is ended with a brief conclusion together with two appendices. 

2 Krylov subspace-Lanczos algorithm 

In this section we give a brief review of some of the main features of Krylov subspace projection 
methods and Lanczos algorithm and more details are referred to [15, 16, 17, 18]. 

Krylov subspace projection methods (KSPM) are probably the most important class of 
projection methods for linear systems and for eigenvalue problems. In KSPM, approximations 
to the desired eigenpairs of an n x n matrix A are extracted from a d- dimensional Krylov 
subspace 



where |c/> ) is often a randomly chosen starting vector called reference state and d <C n. In 
practice, the retrieval of desired spectral information is accomplished by constructing an or- 
thonormal basis Vd G R nxd of Kd(\<f>o),A) and computing eigenvalues and eigenvectors of the 
d by d projected matrix H d = P Vd T AP Vd , where P Vd is projection operator to ci-dimensional 
subspace spanned by the basis Vd- 

The most popular algorithm for finding an orthonormal basis for the Krylov subspace, is 
Lanczos algorithm. The Lanczos algorithm transforms a Hermitian matrix A into a tridiagonal 
form iteratively, i.e., the matrix A will be of tridiagonal form in the ci-dimensional projected 
subspace H d . In fact, the Lanczos algorithm is deeply rooted in the theory of orthogonal 
polynomials, which builds an orthonormal sequence of vectors (|c/>o), \4>i), ■ \<f>d-i)} and satisfy 
the following three-term recursion relations 



K d (\<j> ),A) = spanil^), A\<j> ), • • • , A^lfo)}, 



(2-1) 



A\4>i) = /3 i+1 \(f) i+1 ) + ailfa) + Pi\<f>i-i). 



(2-2) 
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The vectors \<f>i),i = 0, l,...,d — 1 form an orthonormal basis for the Krylov subspace 
K d (\(J) ), A). In these basis, the matrix A is projected to the following symmetric tridiagonal 
matrix: 

/ 

a (3i 

Pi "i Pi 
Lj = p 2 a 3 p 3 



\ 



\ a d -i J 

where the scalars p i+ i and ctj are computed to satisfy two requirements, namely that \4>i+i) be 
orthogonal to |0j) and that |||0i+i)|| = 1. 

In fact, the Lanczos algorithm is a modified version of the classical Gram-Schmidt orthog- 
onalization process. As it can be seen, at its heart is an efficient three-term recursion relation 
which arises because the matrix A is real and symmetric. 

If we define the Krylov matrix K such that the columns of K are Krylov basis {A l (f) ; i = 
0, d — 1} as: 

ir:=(|0 o ),A|0o),...,^ 1 |0o)), 

the application of the orthonormalization process to the Krylov matrix is equivalent to the 
construction of an upper triangular matrix P such that the resulting sequence $ = KP satisfies 
$t$ — i. We denote by \<pj) and Pj respectively the j-th column of <3> and P. Then we have 

(001^(^)^(^)100) = (KPtlKPj) = (0 4 |0,), (2-3) 

where = a + a\A + ... + a^A 1 is a polynomial of degree i in indeterminate A. 

In the remaining part of this section we give an algorithmic outline of the Lanczos algorithm, 
where it will be used in calculation of amplitudes of continuous-time quantum walk. 

Lanczos algorithm 

Input: Matrix A G K nxn , starting vector |0 O ), |||0o)|| = 1, scalar d 
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Output: Orthogonal basis {\</>o), \(fid-i)} of Krylov subspace Kd(\</>o), A) 

/?O = O,|0o) = |0)/|||0)|| 

for 2 = 0,1,2,... 

\v t ) = A\4>i) 

OLi = (4>i\Vi) 

\v i+1 ) = \vi) - Al^i-i) - 
</ 

Pi+i + o 
= \v i+ i)/(3 i+ i 
else 
10m) = 0. 

3 Graphs, adjacency matrix and generalized stratifica- 
tion 

In this section we give a brief outline of some of the main features of graphs such as adjacency 
matrix, distance function and then by turning the graphs into a metric space based on distance 
function, we have been able to generalize the stratification introduced in [13]. A graph is a 
pair T = (V,E), where V is a non-empty set and E is a subset of {(a, (3); a, (3 eV,a^ f3}. 
Elements of V and of E are called vertices and edges, respectively. Two vertices a,/3 e V are 
called adjacent if (a, (3) G E, and in that case we write a ~ f3. For a graph V = (V,E) we 
define the adjacency matrix A by 
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1 if a ~ 




otherwise. 



Obviously, (i) A is symmetric; (ii) an element of A takes a value in {0, 1}; (iii) a diagonal 
element of A vanishes. Conversely, for a non-empty set V, a graph structure is uniquely 
determined by such a matrix indexed by V. 

The degree or valency of a vertex a G V is defined by 



where | . | denotes the cardinality and n(a) is finite for all a G V (local boundedness). A 
finite sequence ao,a±, ...,a n G V is called a walk of length n (or of n steps) if otk-i ~ for 
all k — 1,2, ...,n. For ct 7^ /5 let d{ot,(3) be the length of the shortest walk connecting a and 
p. By definition <9(a, a) = for all a G V and d(a,P) = 1 if and only if ct ~ /3. Therefore, 
graphs become metric space with respect to above defined distance function d. 

Now, in the remaining part of this section we try to define generalized stratification based 
on distance function. To this aim, similar to association scheme [19] we define a partition 
(called distance partition) on V x V, i.e., V x V = \Ji Ti based on distance function <9, where 
the subset are defined by 



Using above distance partition one can define the set T^a) (i-th neighborhood of vertex a) as 



K (a) = \{PeV;a~P}\, 



Ti = {(a,P) G V x V\d(a,P) = i}. 



(3-4) 



T t (a) = {PeV\(a,P)eT l }. 



(3-5) 



Obviously the class of subsets T^a) defined above partition V as 



(3-6) 
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(see Fig.l). As we see the graph is stratified into a disjoint union of strata, hence we call it 
the generalized stratification based on distance function with respect to vertex a, where the 
vertex a is referred to as a reference state(see Fig.2). 

In this stratification for any connected graph T , we have 



for each (5 G T^a). 

Obviously above relations are similar to those of distance regular graphs [19], where in the 
later case the sets form an association scheme and the stratification Ti(a) is independent 
of reference state a, but in an arbitrary graph the generalized stratification depends on the 
choice of reference state. 

In order to study continuous-time quantum walk on a given graph via stratification, we 
define in the following section a Hilbert space which is suitable for Lanczos algorithm. 

4 Hilbert space for the generalized stratification 

Hereafter, we fix a point o e V as a reference state of the graph then with each stratum T k (o) 
we associate a unit vector \<f> k ) in l 2 (V) called unit vector of Axth stratum. The closed subspace 
of l 2 (V) spanned by {|0fc)} is denoted by A(r). In section 6, we will deal with continuous-time 
quantum walk, where A(r) will be referred as walk space denoted by V wa ik, i.e., the strata 
{|(/>fc)} span a closed subspace, where the quantum walk remains on it forever. 
Since {|0/c)} become a complete orthonormal basis of A(r), we often write 



k 

Then for each stratum Y k (o) of generalized stratification, the unit vector associated to /c-th 
stratum is defined by 




(3-7) 




(4-8) 




(4-9) 
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where, \k,a) denotes the eigenket of a-th. vertex at the stratum k and integers g kja > 1 for 
each a G T k (o). We refer to a graph as GQD graph if the coefficients g k:Ct satisfy conditions 
appearing in appendix A (the conditions (A-ii) through (A-iv)). 
By choosing g ktCe = 1 for each a G T k (o), Eq.(4-9) reduces to 

\<t>k) = 7 = N , E M>, ( 4 - 10 ) 



|r fc (o)| aer fc (o) 

where, /c = 0, 1, 2, ... correspond to unit vectors of QD graphs of Ref.[l]. 

In the following we show that, for the QD type graphs the unit vectors of strata given in 
Eq.(4-10), are the same as the orthonormal basis produced via Lanczos algorithm (this is true 
for GQD graphs too, where its proof is referred to appendix A.). To do so, let us consider the 
action of adjacency matrix A over \<p k ) as 

A \^) = ^===: E A %a) 
yJ\T k (o) \ a er k (o) 

1 E E \k + i,u) + -=l= E E IM 



^k(o)\ aer k (o) ver k+1 (o),u~a J\T k (o)\ a er k {o) ver k (o),v~a 



E E l*-i^> 



W°)l 



\f\Fk(o) \ aer k (o) i/ernW^Q M l rfc ( )l ^/l r fc+i(°)l ^er fc+1 (o) 

1 



E A fc+ i(i/)|A; + 



E afc(v)|A;,i/) + 



r*_i(o)| ; E A^(^)|r fc (o)l |A; _ l ^ (4 _ n) 



^1^,(0)1^(0) ' \ ^^(o)! , 6 r fc _ l(o) 

By defining f3 k = jj^pg A fc (i/), A fc (i/) = |{a G r fc _i(o); a ~ and a fe = G v ~ a}| 
for a, z/ G 1^(0), the three-term recursion relations (4-11) reduce to those given in (2-2). 

Therefore, the adjacency matrix takes a tridiagonal form in the basis \<p k ) (orthonormal 
basis associated with strata), consequently these basis are identical with the orthonormal basis 
produced by Lanczos algorithm. 
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5 Spectral distribution of the adjacency matrix A 

It is well known that, for any pair (A, \4>o)) of a matrix A and a vector |0 O ), it can be assigned 
a measure /i as follows 

Ai(a;) = <0o|£(a;)|&>), (5-12) 

where E(x) = J2i \ui)(ui\ is the operator of projection onto the eigenspace of A corresponding 
to eigenvalue x, i.e., 

A = J xE(x)dx. (5-13) 
It is easy to see that, for any polynomial P(A) we have 

P(A) = J P(x)E(x)dx, (5-14) 

where for discrete spectrum the above integrals are replaced by summation. 

Actually the spectral analysis of operators is an important issue in quantum mechan- 
ics, operator theory and mathematical physics [20, 21]. As an example /i(dx) = \ip(x)\ 2 dx 
(fj>(dp) = \ip{p)\ 2 dp) is a spectral distribution which is assigned to the position (momentum) 
operator X(P). Moreover, in general quasi-distributions are the assigned spectral distribu- 
tions of two hermitian non-commuting operators with a prescribed ordering. For example the 
Wigner distribution in phase space is the assigned spectral distribution for two non-commuting 
operators X (shift operator) and P (momentum operator) with Wyle-ordering among them 
[22, 23]. 

Here in this paper we are concerned with spectral distribution of adjacency matrices of 
graphs, since the spectrum of a given graph can be determined by spectral distribution of its 
adjacency matrix A. 

Therefore, using the relations (5-12) and (5-14), the expectation value of powers of adja- 
cency matrix A over starting site |0 O ) can be written as 

(0o|^ m |0o) = / x m u,(dx), m = 0,l,2,.... (5-15) 

JR 
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The existence of a spectral distribution satisfying (5-15) is a consequence of Hamburgers the- 
orem, see e.g., Shohat and Tamarkin [[24], Theorem 1.2]. 

Obviously relation (5-15) implies an isomorphism from the Hilbert space of generalized 
stratification onto the closed linear span of the orthogonal polynomials with respect to the 
measure \i. Since, from the orthogonality of vectors \<f>j) ( Hilbert space of generalized strati- 
fication) produced from Lanczos algorithm process we have, 

6 ij = {<f> i \<t> j ) = (MPi(A)P j (A)\<t>o) 

= J P*(x)P j (xMx)dx = (P^Pj),. (5-16) 

Conversely if Pq, P n -\ is the system of orthonormal polynomials for the measure \i then the 
vectors 

|0,) = P,(A)|0 o ), (5-17) 

will coincide with the sequence of orthonormal vectors produced by the Lanczos algorithm 
applied to (A, \(f>o)). 

Now, substituting (5-17) in (2-2), we get three term recursion relations between polynomials 
Pj(A), which leads to the following three term recursion between polynomials Pj(x) 

f3 k+1 P k+1 (x) = (x- a k )P k (x) - f3 k P k ^(x) (5-18) 

for k — 0, n — 1. 

Multiplying by (3\...(3 k we obtain 

/3i...& + iPfc + i(aO = (x- a k )(3 1 ...(3 k P k (x) - (3 2 k .[3 1 ...[3 k _ 1 P k _ 1 (x). (5-19) 

By rescaling P k as P' k = /3i.../3 k P k , the spectral distribution \i under question is characterized 
by the property of orthonormal polynomials {P^} defined recurrently by 

P&x) = 1, P[{x) = X, 

xP' k {x) = P' k+1 {x) + a k P' k {x) + P 2 k PLi(x), (5-20) 
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for k > 1. 

If such a spectral distribution is unique, the spectral distribution jj, is determined by the 
identity: 

G M= fB^ = ! j, = @M = ±-*!-, (5-21) 

2 " Q 1 "si 



where G^z) is called the Stieltjes transform of spectral distribution [i and polynomials {Q k ^} 
are defined recurrently as 

Q ( o\x) = l, Q?\x)=x-a 1 

xQ ( k\x) = Q k %{x) + a k+1 Qt\x) + Pl^Q^Ux), (5-22) 
for k > 1. The coefficients A t appearing in (5-21) are the same Guass quadrature constants 
which are calculated as 

A l =]im(z-x l )G li (z), (5-23) 

where, Xi are the roots of polynomial P' n (x). 

Now let G^{z) is known, then the spectral distribution /i can be recovered from G^(z) by 

means of the Stieltjes inversion formula as 

1 rv 

fi(y) — fi(x) = lim / Im{GJu + iv )}du. (5-24) 

7T v >0+ Jx 

Substituting the right hand side of (5-21) in (5-24), the spectral distribution can be determined 
in terms of xi, I — 1, 2, ... and Guass quadrature constants Ai, I = 1, 2, ... as 

f i = Y,A l 6(x-xi) (5-25) 
i 

( for more details see Ref. [14, 24, 25, 26]). 

Finally, using the relation (5-17) and the recursion relations (5-20) of polynomial -Pfc(x), 
the other matrix elements (4>k\A m | O ) can be calculated as 

(^\A m | 0o> = * R I x m P^(x)fx(dx), m = 0, 1, 2, .... (5-26) 
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6 Investigation of Continuous-time quantum walk on an 
arbitrary graph via spectral distribution of its adja- 
cency matrix 

Our main goal in this paper is the evaluation of probability amplitudes for continuous-time 
quantum walk by using Eq.(5-26), such that we have 

q k (t) = (<t> k \e~ iAt | O > = * R I e-^P^xMdx), (6-27) 

P1P2 •••Pk JR 

where quit) is the amplitude of observing the walk at stratum k at time t. The conservation 
of probability J2k=o I Qk(t) f = 1 follows immediately from Eq.(6-27), simply by using the 
completeness relation of orthogonal polynomials P' k {x). 

Investigation of continuous-time quantum walk via spectral distribution method, pave the 
way to approximate infinite graphs with finite ones and vice versa, simply via Gauss quadrature 
formula, where in cases of infinite graphs, one can study asymptotic behavior of walk at large 
enough times by using the method of stationary phase approximation (for more details see 

[I])- 

One should note that, the spectral distribution is Fourier transform of the amplitude of 
observing the walk at starting site at time t, i.e., 

q Q (t) = J e- ixt n(x)dx i — ► fi(x) = J e ixt q (t)dt. (6-28) 

Above relations imply that 

q k {t) = / / P' k {x)e- ixt ^x)dx = — 1 / P^x)q (t')e-^dt'dx, (6-29) 

P\P2 ■ ■ ■ PkJ ^P\P2 ■■ ■ Pk J 

therefore, the amplitudes qk(t) can be written in terms of the amplitude qo(t). 
Obviously for finite graphs, the formula (6-27) yields 

Qk(t) = * R E Ae-^P'^), (6-30) 

P\P2 ■ ■ ■ Pk i 
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where by straightforward calculation one can evaluate the average probability for the finite 
graphs as 



In appendix I of Ref. [1] it is proved that for QD graphs the amplitudes on the vertices belonging 
to the one stratum is the same, hence the probability of observing the walk at a site belonging 



appendix I of Ref. [1] is not true any more, consequently the probability amplitudes of observing 
walk at sites can not be obtained from those of strata in a simple way and reader can follow 
the details of calculation of amplitudes in appendix B. 

7 Examples 

7.1 Generalized QD graphs 

Here in this subsection we give examples of GQD type graphs. These graphs look like 
kite and they are embedded in Z k ,k = 2,3,... lattices and defined as follows: Let K(k,n) 
be an k-dimensional lattice graph with n + 1 generalized strata, which consists of vertices 
(0, 0, 0), (I, 0, 0), (0, 1, 0, 0), (0, 0, ...0, /) and (I, I, /) only for odd values of I, where 

k k 

I = 0, 1, n. The vertex (0, 0, 0) is connected to vertices (1, 0, 0), (0, 1, 0), (0, 0, 1) 

, the vertex (0, 0,^/^, 0, 0) is connected to vertices (0, 0, / — 1, 0, 0) and (0, 0, 0, 0) 

i i i 

for each i — 1, k, but for odd values of /, there is an extra connection between (0, 0, I ,0, 0) 



and (/, /, /) (see Fig. 3 for k — 2, n — 6). 

Now, we define unit vectors of generalized strata in such a way that, they coincide with 
the orthonormal basis produced by lanczos algorithm (see appendix A) 




(6-31) 



to stratum k is equal to 



\Qk(t)\ 2 



. Unfortunately for non-QD graphs the lemma appearing in 



|0o) = |O,O,...,O) 



k 
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l<M = 4f: E |1,0,...,0> 

V ft perm. 

\<h) = , 1 A E 12, o, o) + fc|i, i, i)) 



k(k + l) 



perm. 



?2l 



I02i-i) = 4= E |2Z-1,0,..,0) 

V ft perm. 

=L=( 2 |2Z, 0, .., 0) + ft|2Z - 1, 2Z - 1, 21 - 1)), (7-32) 

' ft (ft + 1) perm. 

where, the summations are taken over all possible permutations. Using the relations (A-ii)- 
(A-iv), one can show that the coefficients $ and ctj are 

0[ = k, 0\ = 01 = .... = ft + 1 and a { = 0, i = 1, 2, .... (7-33) 

Now, one can study quantum continuous time walk on these graphs for finite values of n simply 
by following the general prescriptions, but here we restrict ourselves to infinite n. Substituting 
coefficients and ctj in (5-21) , the Stieltjes transform G^z) of spectral distribution u, takes 
the following form 

G,(z) = 1 k • (7-34) 

z 

In order to evaluate above infinite continued fraction, we need first to evaluate the following 
infinite continued fraction defined as 

^ = — An- = — TTT^Tv v t 7 " 35 ) 

z - - I+i - z - (A; + l)G(z) 
z ~ - fc+i 

where by solving above equation, we get 



^ - J* 2 - 4(fc + 1) 

g < 2 » - W) • (7 - 36) 



Inserting (7-36) in (7-35), we get 



G - w = i^kw- (7 - 37) 
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then substituting (7-36) in (7-37), we obtain the following expression for Stieltjes transform of 
H 



(k + 2)z - kJz 2 - 4(k + 1) 

G »« = WTW) • < 7 - 38 ' 

finally by applying Stieltjes inversion formula, we get the absolutely continuous part of spectral 
distribution ji as follows 



= - V \ 2 + l 2 , \x\ < 2VkTT. (7-39) 
Now, we study the probability amplitudes of walk at time t in the limit of large k i.e., 



-iAt 1 f 2Vk+T -ixt k \ ^{k + 1) - X 2 

qi (t) = lim «n\e « I 0o) = lim —== - P^)- V fc2 - ~ 2 dx 



1 / -2 v /(fc+ij7fc . ,/4 (fc + 1 ) k - x 2 

= lim . / , e- lxt P (Vkx)^- 4— dx 7-40 

k ^°° 2n^Jk(k + l) 1 J-^^fik+m " 1 + x 2 /k 

= 1 / 2 e -^ i y (a .)V4^dx = - /* e-^^ 00 (2 3 ;)v / r^rfa;, 
where the polynomial P/^a;) is defined by 



P^x) = lim 1 P/^x). (7-41) 

'-^^(ife + l)' 1 



Now, substituting $ and from (7-33) in three-term recursion relations (5-20), we obtain the 
following relations for polynomials P{(x) 

P^(Vkx) = 1, 

P[(Vkx) = Vkx, 
P^Vkx) = kx 2 - k, 
VkxP[(Vkx) = P{ +1 (\fkx) + (k + ^Pi^Vkx), 1 = 3,4, ... . (7-42) 



Then dividing left and right hand sides of the recursion relations in (7-42) by \Jk{k + 1)'" 
and taking the limit at k — > oo, one can obtain the following recursion relations for P/^rr) 

Pq,oo{ x ) = 1' 
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P-i aolx) = hm — p- = x, 

1 ' ooV ^ fe-oo ^/jfc 

^ (x) = lim -=== = x 2 - 1, 

fe ^°° ^jfe(Jfe + 1) 

^ >00 (x) = J f^ 1>00 (x)+Jt 1>00 (x), « = 3,4, ... . (7-43) 

By comparing the recursion relations (7-43) of P( j00 (x) with those of Tchebichef polynomials 
of second kind, we conclude that 

P^ix) = Ufa/2), (7-44) 

where, C/j(x)'s are Tchebichef polynomials of second kind. Therefore the probability amplitudes 
in Eq.(7-40) can be rewritten as 

2 /-i , 2 fi 



Ql 



(t) = - f e-^U^VT^dx = - [ e~ 2lxt sin((/ + 1) cos" 1 x)dx. (7-45) 

7T J-l 7T J-l 



Now, by doing the change the variable x = cos 9, the integral (7-45) can be written as 

qi (t) = - /V 2l ' cos9 sin((/+ I)9)sm9d9. (7-46) 

7T JO 

Then, using the following integral representation of Bessel polynomials 

JAx) = — r e~ lxcos 9 cos 9d9, (7-47) 

7T JO 

the integral in (7-46) can be written as 

q l (t)=i\.J l (2t) + J l+2 (2t)). (7-48) 

Now, from the recursion relations for Bessel polynomials, i.e., 

21 

J l+1 (x) = -Ji(x) - J,_i(x), (7-49) 
x 

we obtain the following expression for the probability amplitudes of walk in the limit of large 
k 

qi (t) = (l + l)i l *±fV, (7-50) 

where the results are in agreement with the corresponding quantum central limit theorem of 
Ref.[3]. 
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7.2 Non-GQD type graphs 



In this subsection we study an example of non-GQD type graphs, those graphs that do not 
possess three term recursion property. In order to obtain spectral distribution of adjacency 
matrix of a give non-GQD graph, we need to find the basis in which the adjacency matrix has 
tridiagonal form. To this aim we have to choose starting site of walk as a reference state and 
then apply Lanczos algorithm to its adjacency matrix. Then by using spectral distribution, we 
will be able to calculate the amplitudes of walk as will be explained in the following example. 

7.2.1 Walk on finite path graph with second vertex as the starting site of the 



Finite path graph P n = {1,2, ...} is a n- vertex graph with n — 1 edges all on a single open 
path [1]. For this graph, the stratification depends on the choice of starting site of walk. If we 
choose the second vertex as starting site of the walk, as it is shown in Fig. 2, the graph does not 
satisfy a three term recursion relations, i.e., the adjacency matrix has not tridiagonal form. 

Therefore, in order to find the basis in which the adjacency matrix has tridiagonal form, we 
have to apply Lanczos algorithm to the adjacency matrix A of the graph P n , where starting site 
|0 O ) = 1 1 ) is chosen as a reference state. Also, the Lanczos algorithm provides the coefficients 
a and (3 from which the Stieltjes transform G^(z) of //, Eq.(5-21) can be calculated. 

Hence, following the prescription of Lanczos algorithm given in section 2, we get the fol- 
lowing results for P n , which are different for even and odd values of n. 
A. n = 2k 

oti — 0, i — 0, 1, 2k — 1, 



walk 





i = 




^ 2fc -! - 7k- 



B.n = 2k + 1 
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Oi = 0, % = 0, l,...,2k- 1, 
$2i = * — 1) •••) k — 1 

foi-i = i = 1, /c, respectively. 

Substituting the coefficients aij and in (5-20) and (5-22), and using (5-21), we get the 
following closed form of the Stieltjes transform of u. 

G » {z) = IWIW (7 " 51) 

where, C/ n 's are Tchebichef polynomials of second kind. Therefore, the roots x t appearing in 
(5-21) are roots of Tchebichef polynomials of second kind, i.e., xi = 2cos(^-). Also, using 
(5-23) we get the following expression for the coefficients A x 

A = ^-sin 2 (^-). (7-52) 
n+1 y n+V y ' 

Thus, spectral distribution is given by 

" = ^TTg sin2 (^l'^- 2cos (^TT»' < 7 - 53 ' 

Then the probability amplitude of the walk at starting site at time t is 

%{t ) = ^£ sin 2(J^) e -^cosV(n+l) ? (7.54) 
71 ~\~ 1 [—^ 72 H - 1 

again one can calculate the other amplitudes by using Eq.(6-30). 

It should be noticed that, for odd n the Lanczos algorithm produces n — 1 orthonormal 
basis, therefore for calculating the amplitudes on vertices we need to construct an extra vector 
orthogonal to the walk space V wa ik- 

Finally, in the limit of large n, the continuous part of spectral distribution /i(x) is obtained 
as follows 

2 f w 

fj,(x) — — / dy sin 2 (2y) 5 (x — 2cos(y)) 

IT JO 



2 f w sin 2 (2y)5(y — arccos(a;/2)) 



7T JO 



2 sin(y) 
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4 z" 71 



/ dy sin(y) cos 2 (y)S(y 



arccos(a;/2)) 



7T 7o 



= — x 2 \ / A — x 2 , 



-2 < x < 2, 



(7-55) 



therefore, the probability amplitude of the walk at starting site at time t is 



1 f 2 



e~ lxt x 2 V4 



— x 2 dx = 



4Ji(2i) 6J 2 (2t) 



(7-56) 



where above result is obtained by making the change of variable x = cos 8, and using the 
integral representation of Bessel polynomials given in (7-47). Similarly, other amplitudes of 
walk can be calculated by using Eq.(6-27). 

8 Conclusion 

By turning the graphs into a metric space based on distance function, we have been able to gen- 
eralize the stratification and quantum decomposition introduced in [13]. Then the continuous- 
time quantum walk on arbitrary graphs are investigated by spectral distribution method based 
on Krylov subspace-Lanczos algorithm. We have showed that both in QD and GQD graphs, 
the unit vectors of strata are identical with the orthonormal basis produced by Lanczos algo- 
rithm. For more general graphs, we have used the Lanczos algorithm to get a basis in which 
the adjacency matrix has tridiagonal form, where it is necessary for determination of spectral 
distribution of adjacency matrix by using inverse Stieltjes transform. We believe that the in- 
troduced algorithm is a powerful and general tool to investigate the continuous-time quantum 
walk on any arbitrary graph. 



Appendix A 

In this appendix we show that in the case of GQD graphs the unit vectors of strata (i.e., 
Eq.(4-9)), are the same as the orthonormal basis produced via Lanczos algorithm. To do so, 
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let us consider the action of adjacency matrix A over \<pk) as 

A \<l>k) = -i==j= S 9k, a A\k,a) = 1 = 9k, a + 

Y 9k, a Y \ k >v) + -j==z= E 9k, a Y (A-i) 



\fl^a9k,a aer fe (o) i/er fc (o),a~i/ \Jl^a9k,aa&T k {o) i/er M (o),«~i/ 

now in order to have a GQD graph the coefficients g^ a should satisfy the following conditions 



for all v G Tk+i(o) and a ~ v, 



for all i/ G Tfc(o) and a ~ v, 



51 ^fc,a = lk+i9k+i,u, (A-ii) 
aer fc (o) 



X! 9k,a = Vk9k,u, (A-iii) 
aer fc (o) 



2fc,a = 7fc(^ (A-1V) 

for all i/ G Tk-i(o) and a ~ za One should note that the constants 7& and rjk depend only on 
strata number. 



/ J2 9k 

Then by defining f3 k = j k J ^ ° erfc(o) 2 '° and a fc = % for all a G r fe (o),f G r fc _i(o) and 

V ^eer fe _ 1 (o)»fc-i,« 

£ ~ a, the three-term recursion relations (A-i) reduce to those given in (2-2). 

Therefore similar to the QD case, the adjacency matrix takes a tridiagonal form in the 
basis \4>k) (orthonormal basis associated with strata of GQD graphs), consequently these basis 
are identical with the orthonormal basis produced by Lanczos algorithm. 



Appendix B 

Here in this appendix we first prove that in GQD graphs, the ratio of amplitude of a vertex 
in a given stratum to its coefficient appearing in (4-9) is constant, i.e., is independent of 

9k, a 

a G T fc (o). To do so, let us consider the eigenket \(pk) given in (4-9), it is straightforward to 
see that, the eigenket \<pk) together with the following set of states 

!<,) = -=—- — E — M>> / = i,2,...,|r fe (o)|-i, (B-i) 

y^eryo) |^— j2 aG r fc (o) 9k,a 
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form a set of orthonormal basis for a complex space formed by linear span of eigenkets 
belonging to stratum k where u = e |r fc (o)l . 

The above given states are actually orthogonal to all states of walk space (V w ), since 
the eigenket of other stratum do not contain any of \k,a),a G T k (o). Therefore, e~ tAt \<fi ) is 
orthogonal to set of orthogonal vectors for all / = 1, 2, |r fc (o)| — 1; k = 0, 1, d since it 

is a state which remains in V w for all t. Now, substituting (4-9) in (6-27) and {(j>i i\e~ iAt \(j)o) = 
0,1 = 1, 2, |r fe (o)| — 1, we get the following set of equations for amplitudes of vertices 
belonging to stratum k, 

Qk(t) = i 1 = E Qk^qkAt), (B-ii) 
V^ er fc(°) 9k,u aer k (o) 

1 - uu~^ a 

0= I E — ?*,«(*), z = i,2,...,|r fc ( )|-i, (B-hi) 

where qk,a(t) denotes the amplitude of vertex a G 1^(0). To solve equations (B-ii) and (B-iii), 
first we multiply equations (B-iii) by u lv and sum over / = 1, 2, |T fc (o)| — 1, where by using 
the identity ^g^l-i ^-a) = \T k (o)\5 ai/ , we get for v ^ a 

— = TFTTTi zl — n > for a11 <* e r fc(°)- ( B ~ lv ) 

9k,a \ L k{0)\ u ^ a Qk,v 

Above equations imply that gfc '"^ = gfc '^ = for all a,£ G r fe (o) where, is some 
constant independent of vertices of stratum k, and it can be determined by substituting 
qkAt) = B k9k, a in (B-ii) as 

B k = 1 = q k (t). (B-v) 

\jJ2uer k (o) 9k,v 

Therefore, probability amplitude of observing walk at given vertex is proportional to its coef- 
ficient g k>a and it can be written in terms of amplitude of the k-th stratum quit) as 

q k , a (t) = . 9k ' a = q k (t). (B-vi) 



\]Y<veT k {o) 9k,v 
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In QD graphs we have gk, a = 1, for all a G Tk(o), hence vertices belonging to the same 
stratum, have the same amplitude which is in agreement with the result of appendix I of 



In non-GQD type graphs, the coefficients of unit vectors \4>k) do not satisfy the conditions 
(A-ii)-(A-iv) , and we can not obtain vectors orthogonal to V wa \k by the above explained 
prescription of GQD graphs. Therefore, one should use Lanczos algorithm for obtaining n 
independent linear equations, where the amplitudes of vertices of the graph can be deteremined 
by solving them. Let the Krylov subspace generated by the adjacency matrix A and starting 
site |0o) has dimension d, then we will have d unit vectors of strata produced from Lanczos 
algorithm applied to the pair (A, |0 O )) (one should note that the walk space V wa ik is generated 
by applying the Lanczos algorithm to adjacency matrix and starting site of the walk ). In 
the most cases, the dimension of V wa ik is less than the number of vertices (d < n), i.e., the 
Lanczos algorithm applied to the pair (A, \<f>o)), dose not produce the enough basis, therefore 
for obtaining remaining equations we choose new reference states orthogonal to walk space 
V wa ik and then we apply the Lanczos algorithm to the adjacency matrix with new reference 
states, respectively. In the following, we explain the procedure in details for the following 
example. 



where vertex |1) is considered as starting site of the walk. We apply the Lanczos algorithm 
to adjacency matrix A and starting site \<f> ) = |1), where orthonormal basis and coefficients 
ak,Pk produced from Lanczos algorithm are 



Ref.fl]. 



Example 



We consider tree graph of Fig.4, with six vertices and complete orthonormal basis 



{|1),|2),|3),|4),|5),|6)}, 



I0 2 ) 



V2 




1 



(|5) + |6», 




\<h) 



^(|2> + |3> + |4» 

* (_2|2> + |3> + |4», 
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/3i = a/3, /? 2 = & = \AA «i = a 2 = «3 = «4 = 0, (B-vii) 

respectively. One can straightforwardly show that the corresponding Stieltjes transform of \x 
and spectral distribution are 

z 3 - (1 + V2)z/y/3 



z 



4 -(4 + v / 2>7v / 3 + l' 



H = 0.2851952676(<5(a; - 1.662563892) + S(x + 1.662563892)) 

+0.2148047323(5(a; - 0.6014806445) + 5(x + 0.6014806445), (B-viii) 

respectively, which yield the following probability amplitudes of walk at k-th stratum at time 
t, for k = 0,1,2,3 



q (t) = [ e~ ixt fx(dx) = 0.2851952676 cos(1.662563892t) + 0.2148047323 cos(0.6014806445t). 

JR 

qi (t) = -L J R e- ixt P[(x)n(dx) = ±= J RX e-^^dx) 



= — p[0.2851952676sin(1.662563892f) + 0.2148047323 sin(0.6014806445t)], 
iy3 

<&(*) = ^ J R e- ixt Pi(xUdx) = ±= jy - rf) e -**p(dx) 

= 4=[0.2943408772cos(1.662563892t) - .2943408762 cos(0.6014806445t)], 
v2 

q 3 (t) = ^ J R e-^P^{dx) = ^2 J r (x 3 - ^-^x)e-^n{dx) 

= -^[0.102214289 sin(1.662563892t) - 0.2825324240 sin(0.6014806445t)]. (B-ix) 

Obviously, we need two extra equations for obtaining amplitudes on sites of the graph. Ac- 
cording to the above explained prescription we can consider 

l*>> = t^(|5> - |6», (B-x) 

as new reference state (IV'o) £ Km/fc) an d then by applying Lanczos algorithm to the pair 
(A, \ipo)) we obtain 

llM = ^(|3> " |4», (B-xi) 
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which leads to two following extra equations 

(^o|e- iM |0 o > = 0, 

(V>i|e- iM |0 o ) = 0. (B-xii) 

Now, by solving the above six equations, one can obtain amplitudes of continuous-time quan- 
tum walk on vertices of the graph as 

(l\e~ iAt \<h) = <h(t), 

(2|e- lA '|0o) = ^(gi(t)-v / 2g 3 (t)), 
(3|e-^ | O ) = (4|e-^ | O > = ±( qi (t) + ^<fe(*)), 

(5|e-^ | O ) = (6|e-^ | O > = ^fc(f), (B-xih) 
where qo(t), qi(t), qi{t) and 53 (t) have been given in Eq.(B-ix). 
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Figure Captions 

Figure. 1: Shows the stratification with respect to distance function. 
Figure. 2: Shows the finite path graph of P n , where walk starts at vertex 2. 
Figure. 3: Kite graph with k — 2 and n — 6, where walk starts at vertex (0,0). All vertices 
lying on a given vertical dashed line belong to the same statum 
Figure. 4: Shows the tree graph, where walk starts at vertex 1. 
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